program electric
	real	::	E, xmax, t
	integer	::	i, n
	real	:: r, O
	real :: x, y
	real,parameter :: pi = 3.14159265359
	E = 0.1
	r = 1
	do
	do O = 0, 2*pi, 0.2
		x = 2.0*cos(O)*cos(O)/r**3 - 2.0*sin(O)*cos(O)/r**3
		y = (sin(O)**2 + sin(O)*cos(O))/r**3
		print*, x, y, 0.1*x, 0.1*y
	enddo
	r = 0.1*r
	if(r < 0.000001)exit
	enddo
end program
